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Abstract. The limited angular resolution of 7-ray telescopes prevents a direct identification of the majority of sources de- 
tected so far. This is particularly true for the low latitude, probably galactic, ones only 10 % of which has been identified. 
Most counterparts of the identified low-latitude 7-ray sources are Isolated Neutron Stars (INS), both radio-loud and radio- 
quiet (Geminga-like) objects, which are characterised by an extremely high value of the X-ray-to-optical flux ratio fx/ f opt- 
Therefore, the systematic X-ray and optical coverage of low-latitude unidentified 7-ray sources aiming at high fx / /opt sources 
seems one of the most promising ways to spot INS candidate counterparts. Since low latitude sources are heavily affected by 
the interstellar absorption at both X-ray and optical wavelengths, we have focussed on two middle-latitude, probably galactic, 
GRO/EGRET sources: 3EG J0616— 3310 and 3EG J1249— 8330. These two sources, which could belong to a local galactic 
population, have been selected owing to their relatively good positional accuracy, spectral shape and lack of candidate extra- 
galactic radio counterparts. Here we report on X-ray observations of the two 7-ray error boxes performed with XMM-Newton 
and on their optical follow-up carried on with the Wide Field Imager at the ESO/MPG 2.2m telescope. Less than half of the 
~ 300 sources detected by the X-ray coverage have no optical counterparts. Among those, we have selected few interesting 
sources with fx/ j opt > 100, which we consider promising INS candidates. 
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1 . Introduction low-latitudes and up to 1.5° at mid-latitudes) which frustrated 

the search for counterparts at other wavelengths. Moreover, the 
The nature of the Unidentified 7-ray Objects (UGOs) is ,. .. , . .. .. , .. ■ U1 . e 

' 3 J y ' limited 7-ray statistics made it impossible to perform blind 

5h ' one of the main issues of the 7-ray astronomy. The third • .. 1 • ■ a t v a- a txtc 

rrt ■ periodicity analysis aimed at unveiling undiscovered INSs. 

. . , ComptonGammaRay Observatory (CGRO) EGRET catalogue The distribution analysis of the tj GOs in the third EGRET 

< HartmanetalJ1999|) contains 271 high-energy 7-ray sources catalogue (F i gure Q shows that they can be grouped in a t least 

detected at energies above 0.1 GeV. The high-latitude ones four d iff er e nt populations JGehrels et aljEoOfl loTe^lB 

• ' 1 > 10°) are 191, and 93 of them have been identified 2 Q0ll Eool iRomero & et al.ll2004 ~The presumably galactic 



with blazars i.e. featureless flat spectrum radio-loud AGN UGOs (included those observed at | b | > i °) are equally 

and BL Lac objects J von Montigny et alj | l995|) , while 5 of distributed in 3 popu i ationS; with about 45 sources each. The 

the 80 sources at low latitude have been identified with pul- GRp _ T indudes bright and relatively hard sources near the 

sating Isolated Neutron Stars (INSs), both c lassical radio- ^ < 3°). with a concent ration towards the in- 

pulsars and radio-quiet, Geminga-like, objects < Caraveo|2002t ner s piralarm s , Bhattachar v^tliTEool at di stances of a few 

lThompson| | 2004 . In total, 170 EGRET sources, 74 of which kpc feanbach et al.lll99dlRomero et alJll999h . Many of these 



at low latitudes, have no counterpart at lower frequencies and sources m weU correlated with tracers of star formation , which 

remained unidentified. The identification work had been ham- means m at the ir age sh ould be a few million years at most 

pered mainly by the poor localization (about 1 ° in diameter at feomero || 200lt IgTSEooIi . The GRP-II includes variable 

* This work is based on observations obtained with XMM-Newton, sources at \ b \ > 3 ° which ^ distributed in a sort of spheri- 

an ESA science mission with instruments and contributions directly cal hal ° around the Galactic Center with a radial distribution 

funded by ESA Member States and NASA. Optical observations have equivalent to that of globular clusters dGrenienl200lL l2004. 

been obtained with the ESO/MPG 2.2m at La Silla (Chile) under pro- These sources a re softer and significantly more variabl e than 

grams 68.D-0478(A) and 68.D-0478(B). GRP-I sources JTorres et al.ll200 1 alht iNolan et alj Eo031: they 
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Fig. 1. Galactic distribution of the unidentified 7-ray sources 
from the third EGRET catalogue (Hartman et al. 1999): the two 
black circles show the positions of 3EG J0616— 3310 and 3EG 
J1249-8330 



are presumably older (with an age of the order of a few Gyrs) 
and at distances between 2 and 8 kpc. The last group is the 
Local Gamma-Ray Population (LGRP). It is composed by sta- 
ble sources at 3° < |6| < 30°, which are fainter and softer than 
sources located at lower Galactic latitude and have an evident 
asymmetric distribution north of the G ala ctic center and south 
of the anticenter. lGehrels et alJ l koOOh and lGrenieJ OoOO) sug- 
gested that these sources might be associated with the Gould 
Belt, i.e. a 30 Myr old starburst region, 300 pc in radius, com- 
posed by young massive and late type stars, molecular clouds, 
and expanding interstellar medium. Therefore, LGRP sources 
should be young (i.e. a few Myr old) nearby (1 00-400 pc) an d 
low luminosity (~ irj 32-33 erg s _1 ) objects (Grenier 2004). 
But recently, using an improved interstellar emission model, 
ICasandiian et alJ (|2005) have shown that most of these sources 
can correspond to the emission of clumpy dark clouds that sur- 
round all the molecolar clouds of the Gould Belt; therefore their 
existence is not confirmed. Finally, there is an isotropic popu- 
lation of extragalactic origin, which is characterized by a va- 
riety of spectra and variabilities and includes no more than 35 
sources. 

After the release of the third EGRET catalogue, only a few 
additional identifications have been obtained. From the extra- 
galactic point of view, the multiwavelength ap proach has led 
to a b lazar identification of 3EG J20 16+3657 faalpern et alJ 
l2001bl) . 3E G .12006-2321 ftVallace et alJ 120021) and 3EG 
.T2027+3429 iSguera et alJ l2004 . while for 3EG 1621+8203 
( Muk heriee et alJl2002h. 3EG J 1735- 1500 dCombi et all 2003) 
and 3EG J04 16+3650 iS guera et alJl2005t) a radio-galaxy has 
been proposed as counterpart. Within the Galaxy, 3EG J 1824- 
1514 and 3EG J0241+6103 have been associated with the 
two well-known microquasars LS 5039 (McSwain et al. 2004) 
and LS I +61 303 llMassill2004l ICasares et alJl2005t> . while 
ICombi etalJ {2004) suggested the association between the 
microquasar candidate AX J 1639-4642 and the UGO 3EG 



J1639-4702. In two cases a peculiar early-type binary sys- 
tem have been suggested as counterpar t of an UGO, i.e . 
SAX J0635+0533 for 3 EG .10634+052 1 JKaaret et allE oOO) 
and A0535+26 for 3 EG J0542+2610 komero et alJ l200lh . 
Turning now to the INS family, the search for radio or X-ray 
pulsars has provide d likely cand i dates f or 3EG J0222+4253 
(PSR .10218+42 32. iKuiper etalJ |2002h s l. 3EG J1048-5840 
(PSR B 1046- 58. iKasni et alJ ll200ok 3EG J2021+3716 (PSR 
J202 1+3651, Rob erts et alJ ll2002bh and 3EG J2227+6122 
(RX/AX J2229.0+61 14. lHalpern et all d2001 Jah \ Moreover, 
thanks to the P arkes Multi-beam pulsar Survey (PMS) of the 
galactic plane ( Manchester et al. 2001), other promising pul- 
sar candidates have been found fo r 3EG J 1420-6038 and 3EG 
J1837-0606 fo'Amico et alll200ll) and for 3EG J1013-5915 
(Ca milo et alJl200 II) . However, none of the proposed identi- 
fications could be confirmed owing to the lack of contempo- 
rary radio data. Finally, various X-ray and radio studies have 
pointed to a close relationship bet ween the pulsar wind nebu- 
lae (PWN) and the E GRET sources jRoberts et alJ2002alEoOll 
l2005URobertsl2005l) . 

Apart from the cases of Cen A and of the LMC, all 
the EGRET sources identified so far fall either in the INS 
or in the Blazar class. However, also in view of the High 
Energy Stereoscopic System (HESS) results which are un- 
veiling several classes of very high energy 7- ray emitters 
( Aha ronian et al. 2005a bldlctlBrogan et alJl2005l) . it is hard to 
believe that these two classes exhaust all possible 7-ray source 
scenarios. Nevetheless, it has been argued that rotation pow- 
ered p ulsars should dominate the Ga lactic 7-ray source popu- 
lation ( Yadigaroglu & Romanii ri997l) and that many of those 
should be radio-quiet, since the 7-ray beam is broader than 
the radio one. The classic exa mple of a radio-quiet pulsar is 
Geminga dCaraveo et alJll996l) . which offers an elusive tem- 
plate behaviour: prominent in high energy 7-rays, easily de- 
tectable in X-rays but downrig ht faint in optical, w ith spo- 
radic or no radio emission (Bignami & Caraveol ll996l) . In the 
latest years, the idea that galactic UGOs might be associated 
with Geminga-like objects has been strengthened by the cases 
of three UGOs which were associated with radio-quiet INSs: 
3EG J 1835+59 18, the brightest amo ng the unidentified UGOs 
and considered 'the next Gem inga' (Mirabal & Halpern 2001; 
iReimer et alJl200lt lHalpern et alJl2002L and the two sources 
3EG J00 10+7309 and 3EG J2020+4017, positionally coinci- 
dent with the supernova remnants CTA-1 ( Brazier et al. 1998; 
lHalpern et alJ2004l) and 7-Cygni JBrazier et all ^fib -respec- 
tively. 

The energetic of the Geminga-like objects is not suffi- 
cient to account for the very low latitude, presumably more 
distant, GRP-I sources but could account for several mid- 
dle latitude, rather nearby, objects of the LGRP. Based on 
this rationale, we considered the case of two middle latitude 
(see Figure IB UGOs: 3EG J0616-3310 and 3 EG J1249-8330 
1 La P alombar a et af] ll2005l) . Paper I). Both sources are charac- 
terized by a relatively good positional accuracy, with a 95% 
confidence error circle radius #95 ~ 0.6°, by a power law 
photon index T ~ 2.1 and by a steady emission with fluxes 
above 0.1 GeV of (12.6 ± 3. 2) and (19.9 ± 4.4 ) xl0~ 8 
photon cm~ 2 s _1 , respectively (Ha rtman et aljfl999l) . In both 
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Table 1. Main characteristics of the eight EPIC observations of 3EG J0616 3310 (ID = 1^) and 3EG J1249-8330 (ID = 5-8) 



Field 


Rev. 


Date 


R.A. (J2000) 


DEC (J2000) 


Exposure Time (ks) 


AT H 


Detected 


ID 




(UT) 


h 771 S 


O 1 JJ 


PN 


MOS1 


MOS2 


(10 20 cm- 2 ) 


Sources 


1 


346 


2001-10-29T17:04:09 


06 17 47.1 


-32 55 13.9 


6.8 


11.4 


11.5 


2.7 


50 


2 


341 


2001-10-18T23:53:02 


0617 47.1 


-33 25 13.9 


6.7 


12.0 


12.0 


2.5 


37 


3 


346 


2001-10-29T04:27:17 


06 15 24.1 


-33 25 13.9 


2.5 


7.3 


7.7 


2.4 


32 


4 


346 


2001-10-28T23:26:57 


06 15 24.1 


-32 55 13.9 


1.3 


6.3 


7.7 


2.5 


27 


5 


236 


2001-03-23T12:56:43 


12 57 53.1 


-83 15 01.9 


7.0 


11.2 


11.3 


10.2 


38 


6 


236 


2001-03-23T17:54:20 


12 57 53.1 


-83 45 01.9 


8.2 


11.2 


10.9 


8.4 


51 


7 


239 


2001-03-30T03:50:ll 


12 40 13.1 


-83 45 01.9 


0.8 


2.9 


2.4 


9.4 


7 


8 


239 


2001-03-29T22:28:14 


12 40 13.1 


-83 15 01.9 


8.3 


12.7 


12.9 


11.2 


52 



cases, the lack of radio counterparts down to a limit of a 
few tens of mJ v ( Mattox et al. 1 120011: 1 Tornikoski et al. 2002; 
ISowards-Emmerd et alJl2004l) does not support an extragalac- 
tic identification, while an association with the Gould Belt 
appe ars more likely (th eir existence was confirmed also by 
ICasandiian et alJ ( 120051) ). Thus we have applied the multi- 
wavelength approach successfully used for the identification of 
Gem inga as well as for other UGOs associated w ith radio-quiet 
INS dCaraveoll200lHMukheriee & Halpernll2004 . In the mul- 
tiwavelength approach, we start with the X-rays since they are 
the nearest neighbours to 7-rays and can be used to bridge the 
gap between the poor resolution achievable in 7-rays and the 
standards of optical or radio astronomy. First, an X-ray im- 
age of the 7-ray error-box is taken and the detected X-ray 
sources are cross-correlated with optical and radio catalogues, 
either available in archives or produced from ad hoc observa- 
tions. Next, potential INS candidates are singled out amongst 
X-ray sources with an extremely high value of the X-ray-to- 
optical flux ratio /x//opt- This 'top-down' systematic strategy 
is a time consuming exercise since it usually involves several 
observing cycles with different instruments at different facil- 
ities, both space and ground-based. Thus it can be pursued 
for a large number of 7-ray sources only if a semi-automatic 
procedure is set-up. Our work can be viewed as a prepara- 
tory step for the massive identification w ork wh i ch will be 
needed in th e coming years, whe n AGILE jTavani et alj E003) 
and GLAST ( McEn erv et alJ2 004) will deliver hundreds of new 
7-ray sources. The improved angular resolution of these new 
7-ray telescopes will provide arcmin positioning, thus easing 
considerably the identificati on work. To this aim, a statistical 
approach will be necessary jTorres & Reimerl2005l) . 

X-ray observations and data reduction techniques are pre- 
sented in Sj2j while the X-ray source analysis is described in ||3] 
The optical observations are described in while the cross- 
correlations of the X-ray data with the optical and radio ones 
are described in fj5] In ij6]we discuss the source X-ray/optical 
analysis, while the most interesting sources are presented in 
Finally, in fj8]summary and conclusions are outlined. 

2. X-ray observations 

The error boxes of 3EG J06 1 6 - 3 3 1 and 3EG J 1 249 - 8 3 30 are 
circles of ~ 35 arcmin radius, a value com parable to the fiel d 
of view of the XMM-Newton telescopes (Jans en et alJ l2001). 



Furthermore, with an unprecedented collecting area of ~ 2500 
cm 2 @ 1 keV, a good spectral resolution (~ 6 % @ 1 keV) and a 
rather broad energy range (0.1-10 keV), the European Photon 
Imaging Camera (EPIC), is particularly suited to investigate 
faint sources. Therefore, we covered each EGRET error box 
with four ~ 10 ks EPIC observations (see Paper I), correspond- 
ing to ~ 70 % of the 95 % error circle. They did not cover the 
central part of the error circle, which is not a favoured region 
for the true position of the 7-ray source. Therefore, the proba- 
bility that it was actually covered by our observations depends 
only on their geometry and is ~ 70 %. In each observation all 
the three EPIC f ocal pl ane cameras were active: the PN cam- 
era JStruder et alJl200lh was operated in Extended Full Frame 
mode, while the MOS1 and MOS2 cameras llTurneretalJ200lh 
were operated in standard Full Frame mode. In all cases the thin 
filter was used. The details of the eight observations are listed 
in Table [l] where, for each of them, we report the pointing co- 
ordinates, the 'effective' exposure times of the three cameras 
after the soft-proton rejection ( i|2.H . the galactic neutral hy- 
drogen column density in the pointing direction, and the total 
number of detected sources. 

2. 1. Data Processing 

For each pointing we obtained three data sets, one for each 
camera, which were independently processed through the stan- 
dard XMM-Newton Science Analysis System (SAS) v.5.2. In the 
first step, the XMM/SAS tasks emproc and epproc were used to 
linearize the MOS and PN event files, respectively. In the sec- 
ond step, event files were cleaned up for the effects of hot pixels 
and soft proton flares. 

Hot pixels, flickering pixels and bad columns, which pro- 
duce spurious source detections during an observation, are 
largely removed by the on board data processing software while 
the remaining ones are masked by the XMM/SAS. However, we 
found that several spurious events were still present in the pro- 
cessed data. These have been localized using the IRAF task 
cosmicrays and removed using the multipurpose XMM/SAS 
task evselect. We then filtered out time intervals affected by 
high instrument background induced by soft proton flares (en- 
ergies less than a few hundred keV) hitting the detector sur- 
face. To this aim, we accumulated the light curves of the events 
selected close to the CCDs edges and with energies greater 
than 10 keV to avoid contributions from real X-ray source 
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variability. For the PN data we considered only single pixels 
events (PATTERN = 0) with energies 10-13 keV while for the 
MOS data we considered both single and double pixels events 
(PATTERN < 4) with energies 10-12.4 keV in CCDs 2-7. 
Then, we set a count rate threshold for good time intervals at 
0.15 and 0.35 cts s _1 for the MOS and the PN data, respec- 
tively. By selecting only events within the 'good' time intervals 
we finally obtained three 'clean' event lists for each observa- 
tion, whose 'effective' exposure times are reported in Tabled 
As can be seen, the space weather was not favourable during 
the observation of field 7 (which was so much affected by the 
high particle background that we lost up to 80 % of the integra- 
tion time) and, partly, for the observations of fields 3 and 4. An 
example of the processed EPIC/MOS image is shown in Figure 
2, in the case of field 4. 

2.2. Source Detection 

In order to increase the signal to noise (S/N) ratio for the de- 
tected sources and to reach fainter X-ray fluxes, for each ob- 
servation we used the XMM/SAS task merge to merge the clean 
linearized event lists of the 3 cameras. This was possible thanks 
to the excellent relative astrometry between the MOS and the 
PN (~ 1", a value much smaller than the FWHM of the PSF). 
In such a way we obtained a 'total' event list, which was used 
to generate images in 7 different energy ranges. Namely, we 
considered two standard coarse soft/hard bands (0.5-2 and 2- 
10 keV) and a finer energy division (0.3-0.5, 0.5-1, 1-2, 2- 
4.5, 4.5-10 keV). All images were built with a spatial binning 
of 4'.' 3 5 per pixel, roughly matching the physical binning of the 
PN images (4"/pixel) and a factor about four larger than the one 
of the MOS images (l'.'l/pixel). For each energy band, a corre- 
sponding set of exposure maps (i.e. one for each of the three 
detectors) was generated with the XMM/SAS task expmap, to 
account for spatial quantum efficiency (QE), mirror vignetting 
and field of view. Finally, the individual detector maps were 
merged in order to obtain the 'total' exposure map for the rele- 
vant energy range. Both images and exposure maps were used 
as a reference for the source detection, which was performed in 
three steps: 

1 . For each of the selected energy bands, the XMM/SAS task 
eboxdetect was run in local mode to create a preliminary 
source list. Sources were identified by applying the stan- 
dard minimum detection likelihood criterium, i.e., only can- 
didate sources with detection likelihood -In P > 5, where P 
is the probability of a spurious detection due to a Poissonian 
random fluctuation of the background, were validated. 

2. Then, the task esplinemap was run to remove all the vali- 
dated sources from the original image and to create a back- 
ground map by fitting the so called cheesed image with a 
cubic spline. 

3. Finally, for each of the selected energy bands the task 
eboxdetect was run again, but in map mode using as a 
reference the calculated background maps. The likelihood 
values from each individual energy band were then added 
and transformed to equivalent single band detection like- 
lihoods, with a threshold value of 10 applied to accept or 
reject a detected source. 
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Fig. 2. Processed XMM-Newton EPIC/MOS image of field 
4 (see Table ^ Photons selected from three different energy 
bands are colour-coded: red (0.3-1 Kev); green (1-3 KeV); 
blue (3-10 Kev). 

At the end of the detection procedure, for each observation 
we obtained a master source list, which provides the source 
counts, count rates and detection likelihoods for all sources in 
all the detection bands, together with their image position, sky 
position and overall detection likelihood. 

Unfortunately, even using the maximum number of spline 
nodes (20), the fit performed in step 2 (see above) is not suffi- 
ciently flexible to model the local variations of the background. 
Therefore, it was necessary to correct each background map 
pixel by pixel, measuring the counts both in the cheesed image 
and in the backgro und map its e lf app lying the correction algo- 
rithm described in Bal dietalJ J2002). All sources in the mas- 
ter list were checked against the corrected background maps 
and all their parameters calculated again. Then, for all the se- 
lected energy bands, we applied a minimum threshold of 8.5 
on the corrected detection likelihoods —In P, corresponding to 
a probability P(j) = 2 ■ 10~ 4 that the source count number in 
the energy band j originate from a background fluctuation. This 
implies a contamination of at most 1 spurious source per energy 
band. The revised source master list was then filtered to include 
only sources with —In P > 8.5 in at least one of the selected 
energy bands and manually screened to reject residual false de- 
tections. The final master lists contain a total of 146 and 148 
sources for the 3EG J0616-3310 and in the 3EG J1249-8330 
error box, respectively. 

3. X-ray Analysis 

3.1. Global Statistics 

In order to perform a full statistical analysis, we computed the 
number of detected sources in each of the 7 energy ranges, i.e. 
those ones with P(j) < 2 • 10~ 4 in the energy band j. Of 
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course, most of the sources are detected in several bands: in 
Table|2]we report, for each observation of both EGRET fields, 
the number of detected sources in each energy range and its 
percentage over the total (that is the total number of the sources 
which have been detected in at least one energy band). 

We note that almost all sources are detected between 0.5 
and 2 keV, with half of them also detected between 1 and 2 keV, 
while only very few sources are detected at very high or very 
low energies. The number of sources detected in each energy 
band is, in absolute terms, very different across the 8 pointings 
but, taking into account the uneven effective exposure times 
(as it is evident in the case of the observation of field 7, see 
Table 0, the relative number is constant, especially for those 
pointings associated with the same EGRET error box. However, 
we note that below 1 keV the fraction of detected sources is 
indeed lower for the 3EG J1249— 8330 error box than for the 
3EG J0616— 3310 one, probably because of the difference in 
the neutral hydrogen column density, ~ 10 21 cm" 2 and ~ 2.5 • 
10 20 cm -2 , respectively. 



3.2. Count Rate and S/N Distributions 



The histograms of the source count-rate (CR) distribution for 
the two coarse soft (0.5-2 keV) and hard (2-10 keV) energy 
bands are shown in Figure[3] for the single pointings of the two 
EGRET error boxes. Since the number of sources per CR bin 
decreases below the peak CR, we deduce that our sample be- 
comes incom plete for lower C R values (the same approach was 
used bv KickgTaf etalJ ill 99711 ). Therefore, we assume the CR 
peak values as our completeness limit. In all the pointings of 
the 3EG J0616— 3310 error box, the source count-rate distri- 
butions in the low (0.5-2 keV) energy range are pretty similar, 
with peaks at log CR~ -2.6; on the other hand, there are some 
differences among the distributions of 3EG J1249— 8330. The 
marginal differences in the count-rate distributions between 
the pointings of the two EGRET sources are an obvious effect 
of the larger hydrogen column density. On the other hand, the 
count-rate distributions in the high energy range (2-10 keV) 
are significantly different, with peaks at log CR < -3. 

The histograms of the source signal-to-noise (S/N) distri- 
bution in the whole energy range are reported in Figure |4] for 
all the pointings of the two EGRET error boxes. The two dis- 
tributions peak between 4.5 and 5.5, with only a few sources at 
S/N > 10. 

3.3. Source Hardness Ratios Distribution 

Since the count statistics (usually a few tens of photons) of 
the detected sources is too low to produce significant spectra, 
we have performed a qualitative spectral analysis using the CR 



measured in the seven energy ranges defined above to compute 
four different Hardness Ratios (HRs): 

HR1 = [Ci?(0.5-1) - CR(0.3-0.5)]/[Cii(0.5-l) + OR(0.3-0.5)] 

HR2 = [CR(\-2)- CR(0.5-l)]/[CR(l-2) + GR(0.5-1)] 

HR3 = [CR(2-4.5) - CR(0.5-2)]/[CR(2-A.5) + OR(0.5-2)] 

HR4 = [CR(4.5-10) - Ci?(2-4.5)]/[OR(4.5-10) + CR(2-4.5)] 

The histograms of the HRs distributions are shown in 
Figure [5] for the combined pointings of both EGRET error 
boxes. Most sources have HR1 > 0.5, -0.1 < HR2 < +0.4, 
HR3 < -0.5 and HR4 ~ -1, These results, together with those 
of §3.2 based on the count distributions, suggest that, for both 
error boxes, the source population is characterised on average 
by rather soft spectra. 

To obtain a further indication on the sources spectra, we 
compared the measured HRs with the expected ones comp uted 
for two different templa te spectral models JGiacconi et alJ 
2001; Ba rcons et alJ2002l) . namely: a thermal bremsstrahlung, 
with temperatures kT increasing from 0.5 (left) to 5 keV (right), 
and a power-law, with photon indexes T increasing from 1 
(right) to 2.5 (left). In such a way we could identify the spec- 
tral model more appropriate for a given source and tentatively 
assign its most likely spectral parameters. The values of the 
expected HRs are overplotted as vertical bars in Figure [5] As 
can be seen, the distributions are compatible with a rather wide 
range of temperatures and photon indexes, thus suggesting that 
we are indeed sampling different types of sources. This conclu- 
sion is not surprising, since the two areas are at medium galac- 
tic latitude and therefore are expected to contain both galactic 
and extragalactic X-ray sources. 

3.4. Sky Coverage and logN-logS 

In order to compute the source fluxes, we assumed a template 
power-law spectrum with average photon index V = 1.7. For 
each pointing, we estimated the hydrogen column density Nn 
(see T abled from the relation N H = 4.8 x 10 21 x E(B - V) 
cm" 2 (Zombeck 1990) using the average colour excess E(B — 
V) de rived from the dust maps develope d bv ISch legel et all 
dl998l) . Following the procedure used by Bal di et alJ J2002I) . 
count- rate-to-flux convertion factors (CF) were computed for 
the PN and the MOS cameras individually using their updated 
response matrices and then combined with the effective expo- 
sure times of each pointing to derive the total conversion factor 
CF. For each pointing and for each energy range we used the 
exposure maps of each camera and the background map of the 
merged image to compute the flux limit map, which gives, at 
each celestial coordinate, the minimum flux of a so urce in or- 
der to be detected with a probability P = 2 • 10 -4 jBaldi et alJ 
2002). Then, we used the flux limit maps to derive the total 
sky coverage of both EGRET error boxes. These are shown in 
Figure[6] in the two standard soft (0.5-2 keV) and hard (2-10 
keV) energy bands. 

Figure0shows the cumulative logA r -log5 f distributions of 
the sources detected in the two energy ranges 0.5-2 and 2— 
10 keV. For comparison, we overplotted the lower and upper 
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Fig. 3. Histograms of the source count-rate distributions for each EPIC pointing of the error box of 3EG J0616— 3310 {left, Nl- 
N4) and 3EG J1249-8330 {right, N5-N8), in the energy ranges 0.5-2 keV (red line) and 2-10 keV {blue line). For observation 
7 the distribution peaks at higher count rate (logCR = -2.35) due to the shorter effective exposure time. 
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Fig. 4. Histogram of the source S/N distribution for the 3EG J0616-3310 (left, N1-N4) and 3 EG J1249-8330 (right, N5-N8) 
fields. 
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Table 2. X-ray sources detected in each energy range. 
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Range (keV) 


0.5-2 


2-10 


0.3-0.5 


0.5-1 


1-2 


2-4.5 


4.5-10 


Total 


Field ID 


N(%) 


N(%) 


N(%) 


N(%) 


N(%) 


N(%) 


N(%) 


N 


1 


37 (74) 


18 (36) 


11 (22) 


25 (50) 


29 (58) 


18 (36) 


1 (2) 


50 


2 


28 (75.7) 


8 (21.6) 


6(16.2) 


23 (62.2) 


17 (45.9) 


11 (29.7) 


0(0) 


37 


3 


28 (87.5) 


6(18.7) 


7(21.9) 


14 (43.7) 


21 (65.6) 


7(21.9) 


2 (6.2) 


32 


4 


26 (96.3) 


9(33.3) 


4(14.8) 


11 (40.7) 


14 (51.9) 


11 (40.7) 


1 (3.7) 


27 


3EG J0616-3310 


119(81.5) 


41 (28.1) 


28(19.2) 


73 (50) 


81 (55.5) 


47 (32.2) 


4 (2.7) 


146 


5 


32 (84.2) 


8(21.1) 


4(10.5) 


10 (26.3) 


20 (52.6) 


6(15.8) 


1 (2.6) 


38 


6 


37 (72.5) 


13 (25.5) 


6(11.8) 


16(31.4) 


29 (56.9) 


13 (25.5) 


4 (7.8) 


51 


7 


6 (85.7) 


2 (28.6) 


1 (14.3) 


2 (28.6) 


2 (28.6) 


1 (14.3) 


1 (14.3) 


7 


8 


44 (84.6) 


19(36.5) 


3 (5.8) 


15 (28.8) 


26 (50) 


16 (30.8) 


3 (5.8) 


52 


3EG J1249— 8330 


119(80.4) 


42 (28.4) 


14 (9.5) 


43 (29.1) 


77 (52) 


36 (24.3) 


8 (5.4) 


148 
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Fig.5. Histograms of the HR distributions for sources in the 3EG J0616-3310 (left) and 3EG J1249-8330 (right) error boxes. 
Green bars indicate the expected HRs computed for thermal bremsstrahlung spectra with kT = 0.5, 1, 2 and 5 keV (from left to 
right). Magenta bars indicate the expected HRs computed for power-law spectra with T = 1 , 1.5,2 and 2.5 (from right to left). 



limits of the same distributions as measured at high galactic 
latitude, wher e only extr agalactic sources are expected to con- 
tribute jBaldi et al. 2002). In both cases our source density is 
above the upper limit of the high latitude distribution, expe- 
cially in the soft energy band and at low fluxes. This result 
points to a significant excess of galactic sources, whose frac- 
tion is larger at softer energies and lower fluxes. 

4. Optical Observations 

In order to search for the optical counterparts of our X-ray 
sources, we used the Wide Field Imager (WFI) 1 mounted at 
the 2.2m ESO/MPG telescope at La Silla (Chile), where we 
collected data under programme ID 68.D-0478. The WFI is 
a wide field mosaic camera, composed of eight 2048 x 4096 
pixel CCDs, with a scale of O'.'238/pixel and a full field of view 



of 33'7 x 32f7. As it matches very well the diameter of XMM- 
Newton/EPIC field of view (~ 30 arcmin), WFI can provide a 
complete coverage of the targets' area with a minimum num- 
ber of pointings. Observations in the filters [7(877), 5(842), 
V(843), #(844) and 7/(845) 2 were requested in order to max- 
imize the optical spectral coverage and to optimise the objects 
classification in the colour space. To compensate for the in- 
terchip gaps, for each passband the pointings were split in se- 
quences of five dithered exposures with shifts of 35" and 21" 
in RA and DEC, respectively. Observations were performed 
in Service Mode between December 2001 and March 2002. 
Unfortunately, bad weather conditions as well as scheduling 
problems resulted in a highly incomplete and inhomogeneous 
dataset and in a data quality in some cases far from optimal. In 



http://www.ls.eso.org/lasilla/sciops/2p2/E2p2MAVFI/ 



http://www.ls.eso.org/lasilla/sciops/2p2/E2p2MAVFI/filters/ 
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3EG J0616-3310 3EG J1249-B330 




S (erg cm~* s~ 1 ) S (erg cm" 2 s~ 1 ) 

Fig. 6. The total sky coverage of the EPIC observations of the 3EG J0616-33 10 (left) and 3EG J1249— 8330 (right) error boxes, 
in the energy bands 0.5-2 keV (red line) and 2-10 keV (blue line). 
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Fig. 7. Cumulative logiV-logS 1 distribution of the sources detected in the EPIC observations of the 3EG J0616— 3310 (left) and 
3EG J1249— 8330 (right) error boxes, in the energy bands 0.5-2 keV (red squares) and 2-10 keV (blue dots). For comparison, in 
both diagrams we also show the lower and upper limits of the high latitude \0gN-\0gS in the energy ranges 0.5-2 keV (red solid 
lines) and 2-10 keV (blue dashed lines). 



particular, only four of the eight XMM-Newton fields were cov- 
ered by the WFI, i.e. fields 2, 3 and 4 of the 3EG J0616-3310 
error box and field 8 of the 3EG J1249— 8330 error box (see 
Table|3]for a summary of the observations). 



4. 1. Data Reduction 

The data reduction was performed with a pre-release version 
of THELI, a fully automatic pipeline for the reduction of op- 
tical and near-IR imaging data obtained with single- or multi- 
chip cameras. A detailed description of this pipeline and the 
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Table 3. Summary of the optical observations performed by the Wide Field Imager of the 2.2m MPG/ESO telescope at La Silla. 



Date 


Field 


Filter 


Number 


Exposure 


Average 


Average 


dd.mm.yyyy 


ID 


Name 


of Frames 


Time (s) 


Seeing 


Airmass 


06.03.2002 


2 


U 


5 


2500.0 


1.20 


1.15 


06.03.2002 


2 


B 


5 


1500.0 


1.18 


1.28 


10.02.2002 


2 


V 


5 


2000.0 


1.12 


1.30 


05.03.2002 


3 


u 


5 


2500.0 


1.07 


1.14 


10.02.2002 


3 


V 


5 


2000.0 


1.29 


1.16 


05.03.2002 


3 


R 


5 


2000.0 


0.82 


1.05 


08.03.2002 


3 


I 


13 


3250.0 


1.01 


1.18 


05.03.2002 


4 


u 


5 


2500.0 


1.14 


1.30 


12.12.2001 


4 


B 


5 


1500.0 


1.09 


1.05 


12.12.2001 


4 


V 


5 


2000.0 


1.09 


1.01 


12.12.2001 


4 


R 


5 


2000.0 


1.00 


1.11 


10.02.2002 


8 


B 


1 


500.0 


1.69 


1.70 


11.02.2002 


8 


V 


4 


2000.0 


1.64 


1.73 




(to'ie-SO 1 06 h 16"00> 06 h 15"'3O- 06 h I5 m 00- 06 h 14 m 30' 06 h 14"00" 
Right Ascension 

Fig. 8. Composite BVR image (B = blue; V = green; R = red) 
of field 4 taken with WFI at the ESO/MPG 2.2m telescope. 

algorithms used ca n be found in ISchirmer et alJ .20031) and in 
Erbe n et alJ J2005). Here we describe only the steps relevant 
for the the WFI multi-chip detector, such as astrometric and 
photometric calibration, image dedithering and coaddition. All 
pre-processing steps (debiasing, flat fielding, superflatting etc.) 
are similar to those performed with single chip cameras. 

The astrometric solution was performed in several steps. 
First, for each CCD chip pre-processe d image we run 
the S Extractor object detection algorithm .Bertin & Arn outs 
1996) to generate a catalogue which contained high S/N 
(DETECT JV1INAREA = 10, DETECT.THRESH = 10) non 
saturated objects. Objects classified as stars according to the 
SExtractor parameters were then matched, using the crude WFI 
pointing astrometric solution, with the USNO-A2 astrometric 
catalogue which has a nominal intrinsic absolute astrometric 
accuracy of ~ 0'.'25 3 . 



Thereafter, for each chip all the detected objects within 
the overlapping regions of the dithering pattern were cross- 
matched and the results were piped into two \ 2 minimisation 
routines. The first one fits two dimensional third order polyno- 
mials to the optical distortion of each chip, whereas the second 
one determines the relative photometric zeropoint with respect 
to all other chips. The final image coaddition was performed 
chip by chip using as a reference the computed chip distorsion 
maps and the relative chip flux normalization. 

Before performing this step, individual weight maps were 
created for each exposure, based on normalized skyflats. The 
detection of defect pixels was performed on the individual 
chips by means of SExtractor and Eye 4 . Those pixels were then 
set to zero in the weight map. Thereafter, the sky background 
was modelled and subtracted, and the chips were resampled 
and coadded by means of the IRAF ElSdrizzle routine, using 
the weight maps created previously. The registration of the in- 
dividual chips was accurate within ~ 1/1 0th of a pixel (0'.'02), 
so that we did not introduce any artificial decrease of the im- 
age quality. After mosaicing, the final image was associated 
with its corresponding weight map, containing information on 
how long each pixel has been exposed. This is particular use- 
ful since, due to the dithering pattern and the gaps between the 
CCDs, the total exposure time in the coadded image was highly 
inhomogeneous. Using the method described above, we could 
assign to the coadded mosaic a consistent global astrometric 
solution accurate within 0'.'3 and a relative photometric zero- 
point costant within ~ 0.05 mag. 

Figure 8 shows as example a composite WFI image of field 
4 obtained after co-adding the single-passband BVR, reduced 
images. 

4.2. Optical Catalogues 

Source extraction was performed by SExtractor on the final 
co-added images using the previously created weight maps 
to properly account for the varying depths of the mosaics. 
Our parameter combination (DETECT.THRESH = 2.0 and 



3 http://ad.usno.navy.mil/star/star_catsj-ec.html 



4 http://terapix.iap.fr 
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Table 4. Limiting magnitudes (3 a) of our multiband optical-to-IR photometry. UBVRI values are derived either from the WFI 
observations or, when not available, from the GSC-2.3 Bj, F, N equivalent (in italics). JHK values are taken from 2MASS. An 
hyphen indicates that no measurements are available. The last column lists the total number of objects per pointing after matching 
over all the passbands. 



Field ID 


U 


B 


V 


R 


I 


J 


H 


K 


N 


1 




23.0 




22.0 


19.5 


15.8 


15.1 


14.3 


5688 


2 


23.9 


25.1 


24.7 


22.0 


19.5 


15.8 


15.1 


14.3 


20837 


3 


24.5 


23.0 


24.7 


24.6 


23.3 


15.8 


15.1 


14.3 


32202 


4 


23.8 


25.1 


24.7 


24.5 


19.5 


15.8 


15.1 


14.3 


34093 


5 




23.0 




22.0 


19.5 


15.8 


15.1 


14.3 


11329 


6 




23.0 




22.0 


19.5 


15.8 


15.1 


14.3 


7423 


7 




23.0 




22.0 


19.5 


15.8 


15.1 


14.3 


6820 


8 




24.5 


24.3 


22.0 


19.5 


15.8 


15.1 


14.3 


15578 



DETECT JV1IN AREA = 8 pixels) turned out to be the opti- 
mal one as it maximizes the number of detections yielding lit- 
tle contamination by spurious objects, located around saturated 
stars. For uniformity reasons, we applied the same values for 
all catalogues. Even thought the detection threshhold is admit- 
tedly low, we decided to use this value in order to provide can- 
didate optical counterparts for as many X-ray sources as pos- 
sible. The extracted catalogues were overlayed on the images 
and the counterparts were visually inspected. The percentage 
of spurious objects at the images' borders was less than 1%. 

Since a very high photometric accuracy is not critical for 
the purpose of this work, default zeropoints provided by the 
WFI science operation team 5 were used for the photometric 
calibration. The adopted zeropoints (in the Vega magnitude sys- 
tem) are 21.96, 24.53, 24.12, 24.43 and 23.37 for the U, B, V, 
R, and / filters, respectively. 

Since the optical coverage of the two fields performed with 
the WFI is not complete, both in terms of sky and passband 
(see Table 0, we have complemented our dataset using a, yet 
unpublished, version of the Guide Star Catalogue 2 {GSC-2.3) 
which provides photometry in the Bj, F and N passbands, 
overlapping the Johnson's B, R, and I 6 , down to 3 a limiting 
magnitudes of Bj ~ 23, F ~ 22 and N ~ 19.5, with average 
errors of ~ 0.25 (at Bj ~ 20). In addition, to extend our pass- 
band coverage to the near-IR, w e have used the Tw o Micron All 
Sky Survey (2MASS) catalogue llCutri et alJl2003l) . which pro- 
vides photometry in the J (1.25^), H (1.65 fi) and K (2.65 p) 
passbands down to limiting magnitudes of 15.8, 15.1 and 14.3, 
respectively. 

We then cross-correlated all the available catalogues (WFI, 
GSC-2.3 and 2MASS) to produce master optical-to-IR cata- 
logues for each of the eight fields. The four WFI catalogues 
(relative to fields 2, 3, 4 and 8) were cross-correlated in 
turn with both the 2MASS Point Source and Extended Source 
Catalogues through the Gator www interface 7 . Then, the two 
output catalogues were merged and cross-correlated with the 
GSC-2.3. For the remaining fields, we extracted the corre- 
sponding object catalogues from the GSC-2.3 through the 

3 http://www.ls.eso.org/lasilla/sciops/2p2/E2p2MAVFI/zeropoints/ 

6 see http://www.stsci.edu/ for a more detailed description of the 
GSC-2 photometric system 

7 http://irsa.ipac.caltech.edu/applications/Gator/ 



same interface and we cross-correlated them with both the 
2MASS Point Source and Extended Source Catalogues using 
Gator. In all cases, we used a fixed cross-correlation radius of 
1 arcsec which largely accounts for the uncertainty in the abso- 
lute astrometric calibrations of the WFI catalogues (« 0'.'3), 
of the GSC-2.3 (« 0'.'35 8 ) and 2MASS (« Q/.'25 9 ). The fi- 
nal colour coverage for each of the eight fields is summa- 
rized in Table |4] together with the estimated limiting magni- 
tude in each passband and the total number of single objects 
extracted from the optical-to-IR master catalogues. The same 
catalogues are also used for a multi-colour analysis aimed at 
the characterization of the stellar/galactic populations in the 
fields (Hatziminaoglou et al., in preparation). 

5. Cross-correlations 

5. 1. X-rays vs. Optical Catalogues 

Sinc e the absolute coo rdinate accuracy of XMM-Newton is ~ 
5" jKirsch et alJl2004 . i.e. a factor 10 worse than the astro- 
metric accuracy (< 0'.'5) of our optical/IR data, we first tried 
to improve the accuracy of the X-ray coordinates. After over- 
laying the X-ray positions on the Digital Sky Survey images, 
we found few X-ray sources which could be confidently asso- 
ciated with a single bright optical object and we assumed the 
optical positions as the true ones. Then, using the IRAF task 
geomap, we calculated the correction to be applied to the X- 
ray coordinates (value always < 2") and we applied it to all 
the remaining X-ray sources using the IRAF task geoxytran. In 
the following, we decided to use a conservative value of 5" (i.e. 
~ 3 times our astrometric correction) for the cross-correlation 
radius. 

In Table|5]we report, for each XMM-Newton field, the num- 
ber of X-ray sources with and without candidate optical/IR 
counterparts, as well as the total number of candidate coun- 
terparts. It is evident that the results of the cross-correlations 
are significantly different for the different fields, depending on 
the varying limiting magnitude of the optical coverage (see 
Table 0}. For instance, for most X-ray sources of fields 1 and 
5-7 we found no candidate counterpart within 5", owing to 
the limiting flux of GSC-2.3 which is, on average, a factor 6 

8 http://www-gsss.stsci.edu/gsc/gsc2/GSC2home.htm 

9 http://spider.ipac.caltech.edu/staff/hlm/2mass/overv/overv.html 
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Table 5. Results of the cross-correlations between the list of X-ray sources and the optical/IR master catalogues. For all fields, 
the total number of candidate counterparts is larger than those of the X-ray sources because of multiple matches. 



Field ID 


Detected 


X-ray sources 


X-ray sources 


Candidate 


Reliability 




Sources 


with no counterpart 


with counterpart 


Counterparts 


(1-P) 


1 


50 


27 


23 


26 


84% 


2 


37 


7 


30 


46 


76% 


3 


32 


7 


25 


41 


82% 


4 


27 


3 


24 


40 


70 % 


Total 


146 


44 


102 


153 




5 


38 


18 


20 


21 


71 % 


6 


51 


34 


17 


19 


80% 


7 


7 


5 


2 


2 


81 % 


8 


52 


24 


28 


37 


76 % 


Total 


148 


81 


67 


79 





shallower in flux than the WFI catalogues. Indeed, the frac- 
tion of X-ray sources without candidate counterparts decreases 
drastically for fields 2-4. This is particularly true for field 4, 
probably due to its short X-ray effective exposure time (see 
Tabled which results in the detection of only relatively bright 
X-ray sources, with presumably brighter optical counterparts. 
Conversely, the longer X-ray effective exposure time of field 
8 results in the detection of fainter X-ray sources, presumably 
characterized by fainter optical counterparts. 

The cross-correlation between X-ray and optical cata- 
logues is obviously affected by spurious matches. In order 
to estimate it, we computed from the optical/IR master cat- 
alogues the total number of objects within the areas encom- 
passed by each XMM-Newton fields. Then, we used the relation 
P = 1 — e~ 77r p , where r is the X-ray error circle radius (5") 
and /1 is the surface density of the optical objects (per square 
degree), to obtain the chance coincidence probability between 
an X-ray and an optical object (Severgnini et al. 2005). By as- 
suming an X-ray error circle of 5", for each observation we 
estimated the probability P to vary between 16 and 30 % over 
all the eight fields (see Table[5}, which means that, at our limit- 
ing magnitudes, contamination effects cannot be ignored in the 
evaluation of the candidate counterparts. 

5.2. X-rays vs. Radio catalogues 

All X-ray sources were also cross-correlated with radio cata- 
logues, namely: 

- The Parkes Radio Source catalogue (PKSCAT90), 
which includ es 8264 ra dio object at 5 < +27° 
JWright & Otrupcekll990l) 

- The Parke s-MIT-NRAO source catalogue (PMN) for both 
the Southern and the Zenith surveys: the firs t one reports 
23277 radio sources at -87.5° < S < -37° JWri ght et alJ 
119941) : the second one in cludes 2400 sources with —37° < 
6 < -29° JWright et all 19961) 

- The NRAO VIA Sky Survey catalogue (NVSS), which has 
almost 2 million sources at S > —40° JCondon et all 1998) 

We obtained a positional coincidence with a NVSS object 
for just 3 X-ray sources detected in field 1: NVSS 5228 for 



source XMMU J061759. 1-325850, NVSS 5190/MRC 0616- 
329 for source XMMU J061756.6-324735 and NVSS 4464 
for source XMMU J061721.5-3301 10. Of these, only source 
XMMU J061756. 6-324735 has a candidate optical counter- 
part. Moreover, a marginal spatial correlation with NVSS 2588 
was obtained for source XMMU J061546.9-333347 of field 3, 
which has no candidate optical counterpart. For all the other 
X-ray sources no candidate radio counterpart was found within 
the 5" radius error circles. 

6. X-ray/Optical Analysis 

6. 1. The X-ray-to-optical flux ratio Classification 
Scheme 

Observations performed with several X-ray missions have 
shown that different classes of X-ray emitters have different, 
rather narrow ranges of X-rav-to-optic al flux ratios fx/ f op t 
JStocke et alJll99lt iKrautter et alJll999l) . In particular, such a 
ratio is very high (i.e. > 1000) for INSs while it is lower for all 
the other classes of X-ray sources, with no, or small, overlap- 
ping between the respective range of values. Therefore, we can 
use this parameter in order to reject or retain an X-ray source 
as a candidate INS. For the /x//opt values we considered the 
result s of the Hamburg/ROSAT All Sky Survey JZickgraf et alJ 
l2003l) . which provides the typical range of values for each class 
of celestial sources, as a 'reference' classification scheme. We 
considered the same objects classes, i.e. white dwarfs (WD), 
cataclysmic variables (CV), galaxies (G), cluster of galaxies 
(CG) and active galactic nuclei (AGN). Moreover, we consid- 
ered stars as a single class (S). In the Hamburg/RASS the X- 
ray flux is calculated in the ROSAT energy range 0. 1-2.4 keV, 
assuming a 'typical' spectral shape for each class of sources, 
while the optical flux is based on their B magnitude. 

In our case we calculated the source flux in the XMM- 
Newton energy range 0.3-10 keV, and assumed a common 
emission model for all sources (a power-law with photon in- 
dex r = 1.7 and galactic column density). Moreover, as shown 
in the previous section, the passband coverage of the WFI ob- 
servations is incomplete, with only V available in all fields (see 
Table |4j for 90 candidate counterparts (~ 70 % of the total). 
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Table 6. Range of the expected /x//opt f> ux ratios (min/max value in logarithmic units) for the different object classes and EPIC 
fields. The X-ray fluxes refer to the full XMM-Newton energy range (0.3-10 keV), while the optical fluxes refer to the Johnson's 
B (top) and V (middle) bands and to the Bj band (bottom). For each source class the slight differences in the fx/ / op t values are 
due to the different TVh of the various EPIC fields, which affects the count-rate-to-flux conversion factor. 



Field ID 


Stars 


WDs 


CVs 


Galaxies 


Clusters 


AGN 




min/max 


min/max 


min/max 


min/max 


min/max 


min/max 


1 


-3.90/+0.50 


+ 1.00/+2.70 


-1.02/+1.42 


-0.54/+1.83 


+0.38/+2.11 


+0.25/+1.72 


2 


-3.91/+0.49 


+0.99/+2.69 


-1.02/+1.41 


-0.54/+1.83 


+0.37/+2.10 


+0.26/+1.73 


3 


-3.93/+0.47 


+0.95/+2.65 


-1.04/+1.39 


-0.55/+1.82 


+0.38/+2.11 


+0.23/+1.70 


4 


-3.94/+0.45 


+0.89/+2.60 


-1.06/+1.38 


-0.56/+1.81 


+0.38/+2.11 


+0.22/+1.69 


5 


-3.80/+0.60 


+1.11/+2.81 


-0.91/+1.52 


-0.44/+1.94 


+0.48/+2.21 


+0.35/+1.82 


6 


-3.81/+0.59 


+1.10/+2.80 


-0.92/+1.51 


-0.45/+1.92 


+0.46/+2.19 


+0.34/+1.81 


7 


-3.84/+0.56 


+ 1.04/+2.74 


-0.96/+1.48 


-0.47/+1.91 


+0.46/+2.20 


+0.32/+1.78 


8 


-3.78/+0.62 


+1.12/+2.82 


-0.90/+1.54 


-0.42/+1.95 


+0.49/+2.22 


+0.37/+1.84 



Field ID 


Stars 


WDs 


CVs 


Galaxies 


Clusters 


AGN 




min/max 


min/max 


min/max 


min/max 


min/max 


min/max 


1 


-4.33/+0.87 


+0.65/+3.11 


-1.54/+1.63 


-0.89/+2.04 


+0.07/+2.00 


-0.10/+2.25 


2 


-4.33/+0.86 


+0.64/+3.10 


-1.55/+1.62 


-0.89/+2.04 


+0.06/+2.00 


-0.09/+2.26 


3 


-4.35/+0.84 


+0.60/+3.06 


-1. 57/+ 1.60 


-0.90/+2.03 


+0.07/+2.00 


-0.12/+2.23 


4 


-4.37/+0.83 


+0.55/+3.01 


-1.58/+1.59 


-0.90/+2.02 


+0.07/+2.00 


-0.13/+2.22 


5 


-4.22/+0.97 


+0.76/+3.22 


-1.44/+1.73 


-0.78/+2.15 


+0.17/+2.10 


+0.01/+2.36 


6 


-4.24/+0.96 


+0.76/+3.22 


-1.45/+1.72 


-0.80/+2.13 


+0.15/+2.08 


-0.01/+2.34 


7 


-4.27/+0.93 


+0.69/+3.15 


-1.48/+1.69 


-0.81/+2.12 


+0.16/+2.09 


-0.03/+2.32 


8 


-4.21/+0.99 


+0.78/+3.23 


-1. 42/+ 1.75 


-0.77/+2.16 


+0.18/+2.11 


+0.02/+2.37 




Field ID 


Stars 


WDs 


CVs 


Galaxies 


Clusters 


AGN 




min/max 


min/max 


min/max 


min/max 


min/max 


min/max 


1 


-4.19/+0.43 


+0.74/+2.65 


-1.33/+1.31 


-0.81/+1.72 


+0.1 3/+ 1.91 


-0.02/+1.70 


2 


-4.19/+0.43 


+0.73/+2.64 


-1.33/+1.30 


-0.81/+1.72 


+0.12/+1.91 


+0.00/+1.71 


3 


-4.21/+0.41 


+0.69/+2.59 


-1.35/+1.28 


-0.82/+1.71 


+0.13/+1.91 


-0.04/+1.68 


4 


-4.23/+0.39 


+0.64/+2.55 


-1.37/+1.27 


-0. 82/+ 1.71 


+0.13/+1.91 


-0.05/+1.67 


5 


-4.08/+0.54 


+0.85/+2.76 


-1. 22/+ 1.41 


-0.70/+1.83 


+0.23/+2.01 


+0.09/+1.81 


6 


-4.09/+0.53 


+0.85/+2.75 


-1. 24/+ 1.40 


-0.72/+1.81 


+0.21/+2.00 


+0.07/+1.79 


7 


-4.13/+0.49 


+0.78/+2.68 


-1. 27/+ 1.37 


-0.73/+1.80 


+0.21/+2.00 


+0.05/+1.77 


8 


-4.07/+0.55 


+0.87/+2.77 


-1.21/+1.43 


-0.69/+1.84 


+0.24/+2.03 


+0.10/+1.82 



However, for fields 2 and 4 the B band limiting is deeper than 
the V band one probably owing to the large QE of the detec- 
tor. Therefore, to define our classification scheme, we decided 
to compute the measured fx/ f opt ratio for both the B and V 
bands, while, for candidate counterparts with only GSC-2.3 
photometry, the fx/ f opt ratio was computed from the Bj mag- 
nitude. 

Moreover, for each class of sources we had to correct the 
fx/ f opt ratio found in the Hamburg/RASS by accounting for 
the XMM-Newton detection band, the different assumed spec- 
tral models and the different optical filter. To this aim, we have 
devised a procedure which is described in the AppendixlAl The 
computed fx / fopt ranges, relative to the B, V and Bj pass- 
bands, for the different classes of X-ray sources are listed in 
Table [6] The repor ted values are i n rou gh agre ement with the 
resul ts obtained by iRrautter et al J \ 19991) and bv lZickgraf et alJ 
( 2003 ), except for the WD class. In this case, the discrepancy is 
due to the correction applied to the Hamburg-RASS flux ratio in 
order to obtain the corresponding XMM-Newton one. In fact, in 
the first case a very soft blackbody spectrum (T = 50000 K) was 



adopted, while we consider a rather hard (r = 1 .7), less realistic 
power-law spectrum; moreover, the ROSAT energy range (0.1- 
2.4 keV) is very suitable for the WD detection, which is not the 
case for the considered XMM-Newton one (0.3-10 keV). Both 
the items imply that, if the count rate of a WD is high enough to 
be detected by XMM-Newton, we would assign it an unrealistic 
large flux. 

6.2. Evaluation of the X-ray-to-optical flux ratios 

The values of X-ray-to-optical flux ratios fx/ fopt for all X- 
ray sources are shown in Figure [9] for all the available opti- 
cal passbands. When a single X-ray source has more than one 
candidate counterpart, we report the fx/ fopt value computed 
for each candidate counterpart. In the case of X-ray sources 
with no optical counterpart, we estimated the lower limits on 
the fx I fopt ratio from the limiting magnitudes of the different 
fields (see table @J. Thus, we assumed limiting magnitudes V 
= 24.7 for fields 2-4 and V = 24.3 for field 8 while for fields 1 
and 5-7 we assumed a limiting magnitude of Bj = 23. 
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Fig. 9. X-ray-to-optical flux ratios fx/fopt for the X-ray sources detected in the 3EG J0616-3310 (left) and 3 EG J1249-8330 
(right) fields plotted as a function of the X-ray flux. Different colours refer to magnitudes computed in different bands i.e. V 
(green), B,Bj (blue), R,Rf (red). In case of multiple optical counterparts for a single X-ray source, the fx/fopt value of 
each candidate counterpart is shown. Vertical arrows indicate the lower limit on fx/fopt for sources with no candidate optical 
counterpart down to Bj = 23, V = 24.7 (left) and V = 24.3 (right). The dashed lines at log( fx/ f opt) = 2 represent the threshold 
level for a source to be considered as a possible INS candidate. 



The systematic optical identification of all X-ray sources 
listed in our catalogue will be presented in detail in a future pa- 
per (La Palombara et al., in preparation). However, on the basis 
of our X-ray-to-optical flux ratio classification scheme (see 
Table [6j, it is very likely that sources with log (fx/fopt) < 
-1.5 are stars. The maximum value measured for log {fx/fopt) 
is ~ 1.9 for sources with a WFI candidate counterpart (fields 2- 
4 and 8) and ~ 1.2 for those with a candidate GSC counterpart 
(pontings 1 and 5-7). These values are fully compatible with 
those typical of various classes of X-ray sources, especially 
the extragalactic ones (AGN or cluster of galaxies). Therefore, 
if we assume that all the above X-ray sources are indeed inden- 
tified with one of their putative optical counterparts, no matter 
which one, none of them can be considered a likely INS. Of 
course, at this stage we can not exclude a priori the possibility 
that some of the candidate optical counterparts are just unre- 
lated field objects. In this case, the X-ray sources would turn 
out to be unidentified and their corresponding fx/fopt would 
increase, moving nearer to INS values. Overlooked INS can- 
didates will be pinpointed after the systematic X-ray optical 
identification work, now in progress (La Palombara et al., in 
preparation). 

7. Sources with no candidate optical 
counterpart 

As a first step we decided to perform a selection within our 
sample by considering only the 125 X-ray sources with no 
candidate optical counterpart. To pinpoint more robust INS 
candidates we have to select those sources which have the 
highest value of the X-ray-to-optical flux ratio and, possibly, 
are characterized by a soft X-ray spectrum, and thus might 
be associated to Geminga-like neutron stars. According to 



Table [6j all the typical classes of X-ray sources are character- 
ized by log( fx/ f opt) 2, since only the brightest extragalac- 
tic sources can exceed this flux ratio level (due to the problems 
described in |J6] we ignore the case of the WDs). Therefore, we 
have decided to use log(fx/ /opt) = 2 as a threshold value and 
to select the sources whose flux ratio, taking into account also 
the relevant uncertainties, approaches to this value. In this way 
we can reject all the galactic sources, almost all the galaxies 
and clusters of galaxies and most of the AGNs. This selection 
limits our sample to 9 sources. For illustration purposes, Fig. 1101 
shows the positions of these sources overlaied on the WFI V- 
band images. The main characteristics of these sources are 
summarized in Table where, for each source, we list the de- 
tection energy bands, its flux and the corresponding X-ray-to- 
optical ratio lower limit. We also report the more likely spectral 
parameters derived by comparing the source HRs with different 
template spectral models, namely: a thermal bremsstrahlung, 
with temperatures klV = 0.5, 1, 2 and 5 keV; a blackbody, with 
temperatures kT^;, = 0.05, 0.1, 0.2 and 0.5 keV; a power-law, 
with photon index r = 1 — 2.5. 

As it is seen from Table [7] we have singled out 8 INS can- 
didate counterparts associated to 3EG J0616— 3310. While all 
of them are potential candidates, XMMU J06 1807. 6-33 1237 is 
certainly more interesting as it meets both our selection cri- 
teria, since it is characterized both by a relatively high X- 
ray-to-optical flux ratio (> 1.73 ± 0.31) and by a soft emis- 
sion spectrum (it is detected only below 1 keV). Therefore 
this source appears a promising candidates for a Geminga- 
like INS. On the other hand, XMMU J061429.8-333225 and 
XMMU J061450.2-331501, which are both detected in EPIC 
field # 3, stand out as the only sources with an X-ray-to- 
optical flux ratio clearly greater than 2 (> 2.48 ± 0.24 and 
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Table 7. Main characteristics of the X-ray sources with no candidate optical counterpart and log(fx/ /opt) ~ 2. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


OBS 


SRC 


DETECTION BANDS 


kT br 


kTbf, 


r 


Flux 


fx / fovt 






(keV) 


(keV) 


(keV) 




(10~ 14 cgs) 


(log) 


2 


XMMU J061807.6-331237 


0.5-1 


< 5 


0.2 


>i 


2.28 ± 1.23 


1.73 ± 0.31 


3 


XMMU J061429.8-333225 


0.5-2; 0.5-1; 1-2 


> 2 


0.5 


1-2.5 


12.92 ± 3.67 


2.48 ± 0.24 


3 


XMMU J061450.2-331501 


0.5-2; 1-2 


2-5 


> 0.2 


1.5-2 


5.14 ± 2.21 


2.08 ± 0.27 


3 


XMMU J061526.1-331724 


0.5-2; 1-2 


> 2 


> 0.2 


1-2 


3.22 ± 1.12 


1.88 ±0.25 


3 


XMMU J06 1546.9-333347 


1-2 


> 1 


> 0.2 


<2.5 


1.93 ± 1.03 


1.66 ±0.31 


4 


XMMU J061507.9-330026 


0.5-2; 2-10; 1-2; 2-4.5 


>5 


> 0.5 


<1 


3.67 ± 1.18 


1.93 ±0.24 


4 


XMMU J061557.2-324635 


0.5-2; 1-2 


0.5-5 


0.2-0.5 


2-2.5 


2.95 ± 1.47 


1.84 ±0.29 


4 


XMMU J06 1504.5-330533 


0.5-2 


<5 


> 0.2 


>1.5 


2.06 ± 1.41 


1.68 ±0.36 


8 


XMMU J124642.5-832212 


2-10 


>5 


> 0.5 


<1 


3.62 ± 1.69 


1.77 ±0.20 



Key to Table - Col.(l): EPIC field sequential reference number. Col. (2): Source identification. Col. (3): Source Detection Bands, i.e. the energy 
ranges j where the probability P(j) that the source counts are due to a background fluctuation is lower than 2 x 10~ 4 (see text). Col. (4): 
Estimated range of possible kTf, r temperatures for a thermal bremsstrahlung spectrum (i.e. XSPEC wabs bremss model) compatible with the 
measured HRs. Col. (5): Estimated range of possible kTtt temperatures for a thermal blackbody spectrum (i.e. XSPEC wabs bbody model) 
compatible with the measured HRs. Col. (6): Estimated range of photon index V values for a power law spectrum (i.e. XSPEC wabs pow model) 
compatible with the measured HRs. In all cases, we assumed the galactic Nu estimated in the pointing direction (Table0. The wide range of 
values for the spectral parameters kTb r , KTtt and F are due to the low count statistics and to the large error bars on the derived source HRs. Col. 
(7): Estimated source flux in the energy range 0.3-10 keV, assuming a power law spectrum with V = 1.7 and the galactic JVh in the pointing 
direction. Col. (8): Lower limit (in logarithmic units) of the X-ray-to-optical ratio, assuming the X-ray flux of column (6) and the derived 
optical flux upper limits. 



> 2.08 ± 0.27, respectively). However, both sources are de- 
tected up to 2 keV while they are not seen below 0.5 keV. As a 
result, their hardness ratios correspond to a thermal blackbody 
temperature greater than 0.2 keV or to a power-law photon in- 
dex r = 1 — 2.5. These values would point towards younger 
objects than Geminga-like INSs. Indeed, Becker & Triimper 
(1997) have shown that the power-law components of X-ray 
detected INSs have average photon index of ~ 1.9, com- 
patible with those estimated for sources XMMU J061429.8- 
333225 and XMMU J061450.2-331501. Thus, we consider 
XMMU J061807.6-331237, XMMU J061429.8-333225 and 
XMMU J061450.2-331501 our best INS candidate counter- 
parts to 3EG J0616— 3310. Owing to the shallower optical cov- 
erage of the 3 EG J1249-8330 error box (see table gj, we 
could single out only one possible INS candidate counterpart 
(XMMU J124642.5-832212). The X-ray spectral parameters 
of this source are unconstrained but they suggest a rather hard 
thermal spectrum which would not be compatible either with a 
Geminga-like or with a young INSs. However, the X-ray-to- 
optical flux ratio of this source (> 1.77 ± 0.02), similar to the 
one of XMMU J06 1807.6-33 1237, still makes it a possible INS 
candidate. 

8. Summary and conclusions 

Identifying high-energy 7-ray sources is a difficult and time- 
consuming task, owing primarily to the dimension of the er- 
ror boxes that have to be covered at different wavelengths. We 
have developed a semiautomatic procedure, encompassing ad- 
hoc X-ray space observations as well as optical ground based 
ones, aimed at finding high /x//opt candidate counterparts and 
we have applied it to two middle latitude EGRET sources. We 
have mapped the error boxes of 3EG J0616— 3310 and 3EG 
J1249-8330 with eight 10 ksec XMM-Newton pointings and 



we have detected about 300 X-ray sources between 0.3 and 10 
keV, down to flux limits of — 4 x 10" 15 and ~ 2 x 10" 14 
erg cm~ 2 s -1 in the energy ranges 0.5-2 and 2-10 keV, re- 
spectively. Four of the eight XMM-Newton pointings have been 
covered in the optical with the Wide Field Imager (WFI) at the 
2.2m MPG/ESO telescope (La Silla) down to a limiting mag- 
nitude of V ~ 24.5. For the remaining fields, optical coverage 
down to Bj ~ 23 has been provided by the GSC-2.3. For all 
fields, the 2MASS catalogue was also used to extend the avail- 
able colour coverage in the near IR in order not to miss possible 
very reddened counterparts. For about 50 % of the detected X- 
ray sources we found candidate counterparts which yield X- 
ray-to-optical flux ratios comparable with the typical range of 
values of the known X-ray source classes. From the sample 
of the X-ray sources with no identification we have selected 
9 sources (8 associated to the 3EG J0616— 3310 error box and 
1 to the 3 EG J1249-8330 one) characterized by X-ray-to- 
optical flux ratios greater than 100. Although all the selected 
X-ray sources can be considered viable candidate counterparts 
to the two EGRET sources, we are not yet in the position to pro- 
pose robust INS identifications. However, we have singled out 
three interesting sources (XMMU J061807.6-331237, XMMU 
J061429.8-333225 and XMMU J061450.2-331501) which are 
particularly promising counterparts to 3EG J0616— 3310 and 
certainly worth further optical investigations. 

While we shall pursue the study of our candidate counter- 
parts, we note that the need to characterize hundreds of X- 
ray and optical sources, just to discard them, is an unavoidable 
bottleneck, limiting the efficiency of any multiwavelength ap- 
proach. Since the number of serendipitous, unrelated sources is 
proportional to the surface covered, a significant step forward 
will be possible only reducing the dimension of the 7-ray error 
boxes. The next generation high energy gamma ray telescopes, 
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Fig. 10. Finding charts of the X-ray sources listed in Table0 For each source the 5" radius error circle is superimposed on the 
corresponding F-band WFI image. 



AGILE and GLAST, promise to improve the source positioning, 
thus significantly reducing the uncertainty region associated to 
each source. Smaller error boxes can be covered with less X- 
ray and optical pointings, thus reducing both the observing time 
and the number of sources in need of a thorough characteriza- 
tion. Hopefully, this will mark a turning point in the long strag- 
gle towards the identification of the UGOs both as individual 
sources and as a population. 
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Table A.l. Range of the applicable B — V colours for the se- 
lected classes of objects 



CLASS 


(B-V) 


(B-V) 


REFERENCE 




min 


max 




Stars 


-0.2 


+1.8 


Esa (1997) 


WDs 


-0.3 


+1.6 


Fier^eroii^yiL 0^95) 


CVs 


+0.2 


+2.0 


Henden & Honevcutt (1995) 


Galaxies 


+0.2 


+1.6 


Prusniel & Heraudeau 09981 


Clusters 


+ 1.0 


+1.5 


Chieregato (2005) 


AGN 


-0.6 


+1.6 


Veron-Cettv & Veron (2003J 



Appendix A: Correction of the Classification 
Scheme 



According tolchieregatol J2005h . for each class of sources the 
fx/ f opt in the XMM-Newton is obtained from the correspond- 
ing Hamburg/RASS one as follows: 

/x.XMM jx,Hamburg-RASS fx.XMM 



ft 



h 



/x,Hamburg-RASS 



(A.l) 



We considered separately the corrections due to the differ- 
ent energy ranges and to the different spectral models: 



X.XMM 



/o.3- 



10, PL 



jxHamburg-RASS /o. 1-2. 4, MOD 
/o. 3-10, PL c/o.i-2.4,PL 
fo. 1-2. 4, PL c/o.l-2.4,MOD 



For each of the eight XMM-Newton pointings we estimated 
both the flux ratio /o.3-io,Pl//o.i-2.4,pl and the conversion 
factor c/o.i-2.4.pl- In addition, for each class of sources we 
estimated also the specific conversion factor c/o.i-2.4,mod- 

The /i,xmm//v values were obtained from the 
/x,xmm//b ones according to the relation: 



X.XMM 



_ fx,XMM x fB 

fv fs fv 

where 



(A.3) 



f_B = ^Vcga x 10 - {B -V)/2. 5 = j_ 957 x 10 _(fl-V)/3.5 (A 4) 
fv /v,Vega 

For each class of celestial source we considered the range 
of B — V colours which include 99 % of the total sample (see 
TablelAJl. 

For candidate counterparts with only GSC-2.3 photometry, 
we considered the correction factor for the different blue mag- 
nitude: 



fx.XMM _ fx,XMM fs 
fBj fB fB., 

where: 

/b /s,Ve g 



x 10" 



0.822 x 10" 



(A.5) 



(A.6) 



fBj /flj ,Vega 

and B - Bj = 0.28 X (B — V) dReid et al.lll99lT) . Also 



in this case we referred to Table lA.ll for the B — V colours of 
each class of source. 



Appendix B: Measurement of the 
X-ray-to-optical flux ratios 

We computed the fx.XMM If opt ratios for all the candidate 



counterparts. X-ray fluxes are the same determined in A3 .41 
while optical fluxes were obtained from the measured WFI 
magnitudes using the Pogson formula f m = / m= o x 10~ m / 2 ' 5 , 
assuming f m=0 = 2.904, 6.478, 3.240, 3.828 and 1.997 xlO 6 
erg c m" 2 s" 1 for m = U, B, V, R and I passbands, respec- 
tively llZombeckll990l) . 

Since for the Bj band, there is no simple relation to obtain 
the optical flux, we first estimated its Vega band flux. This is 
given by the relation f Bj =o = fx a{f (Bj) x W cS (Bj), where 
^cS (Bj) and Wes(Bj) are, respectively, the effective wave- 
length and width of the filter, while f\ ca (Bj) is me specific flux 
calibration at A e ff (Bj) for Bj = 0. By referring to the Asi ago 
Database on Photometric Systems jMoro & MunarilEoOOl) we 
found A cff (Bj) = 4731 A and W eS (Bj) = 1333 A. In or- 
der to obtain the specific flux calibration f^ rm , we consid ered 
the spectral distribution of Vega provided by Hayes ( 1985) and 
performed a linear interpolation between 4725 and 4750 A. 
In such a way we obtained /4731 = 5.457 x 10~ 9 erg cm -2 
s^A- 1 and f Bj=Q = 1333 x 5.457 x 10~ 9 = 7.274 x 10~ 6 
erg cm" 2 s _1 . Finally, we used the Pogson formula fs, = 
/b 7= o x lQ~ Bj ' 2 ' 5 to obtain the optical flux for a given mag- 
nitude Bj. In the case of the F band magnitude, we applied 
the same procedure but assuming A c g(F) = 6555 A, W e g(F) 
= 767 A and / 6555 = 1.731 x 10" 9 erg cm 



2 s 1 A thus 



(A.2) obtaining f F = 1.328 x 10" 



10 -f/2.5 ergcm -2 s -i_ 



